Bose-Fermi Mixtures in One Dimension 
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We analyze the phase stability and the response of a mixture of bosons and spin-polarized fermions 
in one dimension (ID). Unlike in 3D, phase separation happens for low fermion densities. The 
dynamics of the mixture at low energy is independent of the spin-statistics of the components, and 
zero-sound-like modes exist that are essentially undamped. 
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Binary mixtures of dilute quantum gases are a subject 
of steadily growing interest initiated by the realization 
of Bose Einstein condensation (BEC) of alkali atoms 
and motivated by the quest for and subsequent exper- 
imental realization of degenerate Fermi gas Strong 
s-wave interactions which facilitate evaporative cooling 
of bosons are absent among spin polarized fermions due 
to the exclusion principle; so the method of choice for 
cooling fermions to degeneracy has been through the me- 
diation of fermions in another spin state Q or via a buffer 
gas of bosons § §. Dc generacy in dilute gases can be 
understood better than in their liquid helium counter- 
parts due to the weaker interactions, and thus they offer 
prospects of detailed quantitative study of some of the 
most interesting phenomena in the physics of many body 
quantum systems such as the Bardeen-Cooper-Schrieffer 
(BCS) transition g. 

On another front a new generation of BEC experiments 
on surface micro-traps Q and experiments on creating 
atomic waveguides have generated interest in quan- 
tum gases in lower dimensions. Effective one and two 
dimensional BECs have been created, in which excita- 
tions in the confined directions are energetically not al- 
lowed Bose condensation on optical lattices || are 
being actively studied by several groups; the atoms at 
each lattice site can be in regimes of effective ID. 

It is therefore a natural step to bring these two excit- 
ing developments together and consider binary mixtures 
of quantum gases in effective ID, with the possibility 
of forming one dimensional degenerate Fermi gases and 
fermionic waveguides. Fermions in one dimension have 
been the subject and the source of some seminal models 
in many-body quantum physics |lO] mainly because they 
are theoretically more tractable than in 3D. Now there 
is actually the possibility of testing some of these mod- 
els experimentally. Considerable recent theoretical work 
has already been done on three dimensional Bose-Fermi 
mixtures fy], |l2[ |lf, [m[ [jj, |J, |o) but very 
little has been said about one dimensional systems. The 
goal of this paper is to study theoretically some of the 
relevant properties of binary mixtures of bosons and spin 
polarized fermions in an effective one dimensional config- 



uration. In particular we will consider their miscibility 
properties, phase stability and their excitations. 

Model: We consider a longitudinally homogeneous one 
dimensional mixture of iVj, hard-core bosons with mass 
nib and Nf spin-polarized fermions with mass mj at 
T = K. A natural choice of trap-geometry to consider 
such a mixture is in a toroidal trap |2l], ^2) with no ex- 
ternal potential along the circumference (of length L), 
but with tight cylindrically symmetric harmonic confine- 
ment of frequency loq in the transverse direction. This 
geometry can equally well be interpreted as an infinitely 
long, straight waveguide with periodic boundary condi- 
tions. For the atoms to have effective ID behavior at 
zero temperature the ground state energy of the trans- 
verse trapping potential has to be much higher than 
the ground state energy of the bosons and the fermions 
in 3D, i.e. huj > Hb(3D) and hoj > e /(3£>) where 
/ i 6(3£>) — 4:irh 2 ab/mb is the bosonic chemical potential 
and £f(3D) — ^ 2 (67r 2 «/(3D)) 2 / 3 /(2m/) the Fermi energy 
of non-interacting fermions, both in 3D. A measure of 
the transverse spatial extent of the atoms in the ring is 
given by the single-particle ground state widths for the 
transverse trap, = yjfi/mk^o where k — > b, f and 
bf correspond to bosons, fermions and twice the reduced 
mass of a boson and a fermion m&/ = 2m frrib / '(m / + nib) ■ 
A torus of high aspect ratio would have L » r^. 

Our treatment of the ID Bose-Fermi mixture will rely 
on an effective Hamiltonian describing the longitudinal 
behavior of the gas in the toroidal trap 



H = dxip 



dx ip 



2m 



- d x - 
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ipf + 9bf / dx ipl^bipjipf 
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Here ipb(x) and tpf( x ) are field operators for the longi- 
tudinal degree of freedom and x is the circumferential 
spatial coordinate. We have assumed factorization of the 
transverse degrees of freedom; such a factorization is jus- 
tified in regimes of effective ID 0] since the transverse 
spatial dependency is that of the single particle ground 
states 4>bo{f) and (jjfQ^r) for the trapping potential re- 
gardless of the longitudinal behavior or statistics. Thus 
our effective ID coupling strengths for the boson-boson 



and boson-fermion g b f interactions are 



/27rr<ir|0£ ) oM| 4 = 2TiojQa b 

,(r)\ 2 \0 fo {r)\ 2 =2hLu a bf (2) 



9b 

m b 
4irh 2 a bf f 
m bf J 



with a b and a b f being the respective scattering lengths. 
The linear density operators are p b (x) = ipl(x)ip b (x) and 
Pf(x) = ^j(x)tpf(x), with spatially constant equilibrium 
expectations n b — N b /L nj — Nf/L; the density fluc- 
tuation operators are therefore 8p b (x) = p b (x) — n b and 
Sp f (x) = p f (x) -n f . 

Phase stability in the static limit: We first consider the 
mixture in static equilibrium in which case the expecta- 
tions of the fluctuation operators are zero; the kinetic 
energy of the bosons vanishes while the kinetic energy 
for the fermions contributes the Fermi energy per particle 
€f = gfn 2 /3, with gf = h' 7r 2 / '(2m f) and the Fermi wave- 
vector kf = nut. In this static case the total number of 
particles is fixed, so we take the ground state expectation 
of the canonical Hamiltonian for the system and thus ob- 
tain a simple expression for the total energy of a uniform 
mixture of bosons and fermions at equilibrium 



E„ —L 



9f 3 

2 n b + Y n f + 9bfn b n f 



9b 2 



(3) 



The first derivative of this with respect to the densities 
yield the Thomas-Fermi equations 



Pb = gbn b + g b fUf 



Pf = gb/m + gfn 2 f , 



(4) 



for the chemical potentials, and the derivative with re- 
spect to the linear- volume (L) gives the pressure 



P 



dE 
~dL 



9b 2 , 
= ~2 Ub 



2gr r, 

—n f +g bf n b n f . 



(5) 



The second derivative condition for a stable minimum 
with respect to small changes in the densities puts a lower 
limit on the fcrmion density 



n f > 



9 2 bf . 

2g f g b ir 2 a b r 2 



l bf 



(6) 



/ 



This constraint is the opposite of that in 3D where the 
stability condition puts an upper limit on the fermion 
density. The reason for the difference is that the power 
law of the density dependency of the Fermi energy 
changes with dimensionality. The energy contribution 
from the Fermi pressure grows faster as a function of lin- 
ear density in ID than it does with increase in bulk den- 
sity in 3D, however the boson-fermion interaction energy 
behaves similarly in 3D and in ID with respect to bulk 
density and linear density respectively; thus at higher 
fermion densities, the total energy in ID is more likely 
to be lowered if the fermions are spread out over a larger 
volume mixed in with the bosons. 
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FIG. 1: Phase diagram for a mixture of bosons and fermions 
in one dimension. The thin line corresponds to the linear 
stability condition in Eq. (Q). 

The stability criterion in (|^) applies for small fluctu- 
ations, we now analyze the general phase stability for a 
Bose-Fermi mixture in ID; such an analysis was done for 
mixtures in 3D by Viverit et al. p3| . A binary mix- 
ture can have at most two distinct phases that we label 
i = 1,2. The linear volume occupied by each phase is Z,, 
the number of bosons (fermions) therein iVf,(n j and the 
corresponding densities = Nuf^i/Li. The volume 

fractions of the phases are I = L\jL and 1 — £ = L2/ L 
and ratio of the densities in the two phases are labelled 
%(/) = n Hf),i/ n b(f)-2- The total energy for the phase- 
separated mixture is 



E s — LiE,, 



E t \9b 2 , 9f 3 , 
L i { 2" n b,i+ y n fA 9bfn b ,ijif,i 



(7) 



Equilibrium between the phases require 

(a) pi=p 2 ; (b) M6(/),i = Mh(/),2 

(c) Pb(f),i > Pb(f),j if nb{f),i = 0, (8) 

where the pressure and chemical potentials in each phase 
are given by Eqs. (^) and (||) with the total densities 
replaced by partial densities. We will use the identity 
^ n b(f),i + (1 — ^)«-&(/),2 = n b(f) in Eq. (||) to evaluate 
the energy E u of the uniform phase to compare with the 
energy E s of the phase-separated mixture in Eq. (0) . It is 
convenient to introduce density measures in terms of the 
interaction strengths: C/ = 9 b f I (9 f9b) for fermions and 
Cfc = 9bf/(9f9b) f° r bosons. There four possible ways of 
phase separation, the feasibility of each is determined by 
the specific nature of the conditions (||) and the principle 
of minimum energy. We now discuss each case, leaving 
out the somewhat tedious algebra for brevity. 

(i) Two pure phases : The fermions are all in one phase 
and the bosons in the other, so we set nt\ — n b 2 = 0. 
The equilibrium conditions ((8^) and (||c) constrain the 
partial densities: n/ t 2 < 3C//4 and n b: i < 3Cf,/4. When 
those conditions are used in Eqs. (||) and (0) they give 



[E u - E s ] > (1 - i)t 



i9f 3 
n F2 



>0, 



(9) 
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which means that the separated phase has lower energy 
for all values of the volume fraction £ £ [0, 1] and hence 
is energetically preferred in the density regimes where 
phase equilibrium is possible: 



n s <{\-£)-C f 



n b < £-C b . 



(10) 



(ii) A mixed phase and purely bosonic phase : The 
fermions are all in one phase, = but there are 
bosons in both phases. The two equations arising from 
pressure equality (0a) and the equality of boson chemi- 
cal potentials (@b) fix the fermion partial density n/,2 = 
3C//4 which obeys the condition n/,2 < Cf due to the 
inequality of the fermion chemical potentials > p,f,2- 
Then it follows from fib.i = pb,2 that n^i — n,b,2 = 3C&/4. 

On applying the equations for the partial pressures and 
the boson chemical potentials to Eqs. (||) and (Q) we find 



[E u - E s ] = (1 - £)(■ 



:9f 



l F,2 



> 0. 



(11) 



so that the separated phase is energetically preferred in 
this case as well, in the density regimes given by 



3C 



nb = n b , 2 + £ 



(12) 



that do not overlap with those of the previous case. 

(iii) A mixed phase and purely fermionic one : All the 
bosons are in one phase, n b _2 = while the fermions 
can be in both phases. The phase equilibrium conditions 
require that the fermion density ratio r/f £ (0, 1) and that 
it satisfies the equation 



3n^ 
2C f 



(1 + Vf) 2 = (2 + »?/). 



(13) 



This limits the fermion density in the second phase to be 
Cf/2 < n/,2 < 4C/-/3. However on using the equality of 
partial pressures (|8ja) and that of the fermion chemical 
potentials in the two phases, — fj,f t 2 we find that 



i[E u -E s ]=£(l-£) 2 ^-[n FA 



n F ,2? < 0, (14) 



where the inequality holds because r]t < 1. The uniform 
mixture will thus be energetically preferred for all values 
of £, hence this phase separation will not occur. 

(iv) Two mixed phases : Finally we consider the case 
where both phases have fermions as well as bosons. There 
are three equations arising from the equilibrium condi- 
tions (|^a) and (||b) which lead to the following equation 
for the fermion density ratio 



(1- Vf f = 



(15) 



with solution rjf = 1 which, due to (ib,i — Pb, 2, also 
implies that the boson density ratio rjb = 1 which means 
that the only allowed solution is when the entire system 
is uniform, and there is no phase separation of this type. 



The allowed phases for different linear densities are 
plotted in Fig. (jl]); phase separation occurs for low 
fermion densities, in qualitative agreement with the lin- 
ear stability condition. Taking the bulk density to be 
ny( 3£) -) ~ nf/irr 2 , the criteria for single phase nj > \Cf 
and effective one-dimensionality £/(3_d) ^ hio® give the 
limits of fermion density for which bosons and fermions 
in ID can coexist in a single phase: a 2 j / \abr 2 j) < rif < 
^/rf. Transverse trap widths r/ ~ 1 /mi achievable 
currently would allow single phase mixtures for fermion 



bulk densities up to n 



10 18 m 3 , which corre- 



sponds to the density of the coldest 6 Li samples in a 
recent experiment |24|] that created a degenerate system 
of bosons ( 23 Na) and fermions ( 6 Li) in 3D. Scattering 
lengths for alkali atoms are of the order of ~ 1 — 10 
nm, that would allow a range of few orders of magnitude 
of fermion density where a stable uniform ID mixture 
of bosons and fermions would form; this range can be 
widened by increasing the transverse trap strength or re- 
ducing the boson-fermion scattering length. Unlike in 
3D, phase separation effects in ID can be observed by 
reducing the density which is usually easier to do than 
increasing it. 

Dynamic response: We now consider the dynamical 
properties of the mixture. For weak interaction strengths 
and low energy modes we can use linear response the- 
ory which is a convenient formulation of first order time- 
dependent perturbation theory in the interaction picture. 
We consider small density fluctuations of the bosons and 
the fermions about equilibrium. The boson fluctuation 
can be considered a density dependent perturbation for 
the fermions and vice versa, so that we have two cou- 
pled linear equations for the expectation of the density 
fluctuations 8pb{x) and 5pf{x) 



8pb(q,u) 



Xb- 9bfSp f {q,uj) 
Xf ■ 9bfSp b {q,uj) 



with retarded density-density response functions 



4ek 



n\6 P '(k)|0}| 2 



n^0 



2lu 



»0 



n0 



(16) 



(17) 



The small imaginary shift ir\ preserves causality and the 
ground state |0) represents the Fermi sea for the fermions 
and the condensate for bosons. In the Bogoliubov ap- 
proximation for the bosons the response function as well 
as the quasi-particle spectrum in ID have algebraic forms 
identical to those in 3D and are given by 



Ub(q) 2 



n b q 



m b [uj 2 -u)%(q)] 
{e q /hf + (v b q) 2 



(18) 



with free quasiparticle energy e q — h q /(2rrib) and 
sound velocity Vb — y 1 g b n bl m b -The poles correspond to 
the energies of the collective modes which are undamped 
in the Bogoliubov approximation. 
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The fermions being spin-polarized do not have s-wave 
interaction so that the response function for the fermions 
is taken to be that for free fermions; unlike the bosonic 
response function this has a form in ID quite distinct 
from that in 3D: 



%2 2 



2Trh 2 q 



In 



(lu + irj) 2 - uj 2 _ 



[to - 



4m 2 



{{k f ±qf-k)f 



(19) 



The calculation leading from the general expression 
Eq. @ to Eq. (|l|) is analogous to that in 3D Q 
with the Fermi-sphere replaced by a "Fermi-interval" 
[— kF,kp]- It is apparent that Im x/(g,w) ^ only if 
< \lu\ < |w+|. 

The low energy, long wavelength collective modes are 
of particular experimental interest, so we do a Taylor ex- 
pansion of the expression ( |l9| ) for the fermionic response 
function Xf for small values of q, but keeping the ratio of 
the energy transfer to momentum transfer uj/q constant. 
The result is quite interesting 



m/[(w + ir]) 2 -ujf(q) 2 } 
u f (q) 2 = (e q /h) 2 + (v f q) 2 . 



(20) 



It is apparent that if we replace the Bogoliubov sound ve- 
locity with the Fermi velocity Vb — > Vf = hkf /m, the real 
part of this limiting form is identical to the Bogoliubov 
density-density response function; and the spectrum cor- 
responding to its poles are identical in form with that of 
the Bogoliubov poles. 

This equivalence of the bosonic and fermionic density 
fluctuations is distinctly a property of one dimension with 
no analog in higher dimensions. Such an equivalence 
is not surprising when one recalls the Luttinger liquid 
model of Haldane |2(| where the low energy behavior of 
quantum fluids in one dimension were shown to be in- 
dependent of spin-statistics. However it is important to 
note that the fermionic response function that we con- 
sider above is that for free fermions while the Luttinger- 
Tomonaga ]lC| model assumes long range interactions 
among the fermions. 

The absence of interaction among fermions distin- 
guishes the qualitative nature of the fermionic excitations 
from the Bogoliubov modes, despite the similarity in the 
low energy structure of the response functions. Strictly 
speaking the fermionic excitations are elementary excita- 
tions, whereas the Bogoliubov modes are collective modes 
of the bosons for which the interactions play a central 
role. One could of course interpret the fermionic exci- 
tations as zero sound modes in the limit of vanishing 
interaction; in that limit we recall that the zero sound 
velocity coincides with the Fermi velocity even in 3D. 

The similarity of the response functions at low ener- 
gies combined with the fact that homogeneous Bogoli- 
ubov modes have the same form in ID and 3D allows us 



to directly apply the results obtained for spatially uni- 
form binary mixtures of bosons in three dimensions to 
the Bose-Fermi mixtures in ID. As an example, we con- 
sider the normal modes of the mixture determined by the 
vanishing of the coefficient determinant of the response 
equations (|l6|), 1 — g 2 fXbXf — 0; which leads to an ex- 
pression for the normal mode velocities similar to that 
for binary mixtures of bosons in 3D 



0* +Vf)± J(vb-v f ) 2 +- 



n f n b 

hf 
1 mfrrib 



(21) 



Here we have used the linear dispersions Wb{q) — Vbq and 
u)f(q) ~ vjq for long wavelength modes. Hydrodynamic 
equations for boson-fermion density fluctuations in the 
collisional regime also give a similar expression for the 
sound velocities [B , but with the crucial difference that 
the fcrmion sound velocity in that case is that of first 
sound. In the static limit q — > we find that the condi- 
tion v± > 0, necessary for positive compressibility, leads 
to the same condition obtained earlier in Eq. (|6[) from 
energy considerations. 

The frequencies corresponding to the original low 
energy Bogoliubov phonons shift due to the interac- 
tion with the fermions by about Slo = uj q — Vbq — 
n f n b9bf1 2 /mbTn f(v 2 — v 2 ); the shift is positive or neg- 
ative depending on whether Vb > Vf or Vb < Vf, simi- 
lar to the behavior in 3D Jy|. But these modes in the 
mixed system differ from their analog in 3D mixtures 
in that they are not damped unless the velocity matches 
the Fermi velocity as seen from Eq. (p0[). In 3D mixtures, 
such modes are damped if Vb < Vf . 

For modes with higher momentum where the exact 
fermion response function ( |l9| ) has to be used, there is 
damping for mode frequencies in the range \u)—\ < M < 
with the damping rate given to lowest order by 
7 ~ mcbg 2 f/(4h 2 gb). The key difference from 3D is that 
this rate is independent of the mode. 

The exact fermionic response function (^) in one di- 
mension has several interesting features distinct from 3D. 
For zero energy transfer it is seen that X/(<Z) 0) nas a 1°S _ 
arithmic divergence at q = 2k f due to perfect nesting, 
whereas in 3D the derivative of the response function 
is divergent, which leads to Friedel oscillations. In ID 
Bose-Fermi mixtures, the logarithmic divergence of re- 
sponse function leads to periodic density variations in 
the fermions of period 2k f associated with the forma- 
tion of coherent superposition of particle- hole pair states 
called the Peierls channel |^7| . Due to the boson-fermion 
density coupling, the bosons acquire a similar periodicity 
but out of phase with the fermion density modulation as 
demonstrated variationally by Miyakawa et al. p8| . 

In conclusion, we have studied the phase stability of 
a boson-fermion mixture in one dimension and demon- 
strated that phase separation would occur at low fermion 
densities a behavior opposite to that in 3D. This means 
that phase separation effects may be studied at densi- 
ties easier to achieve than in 3D. The regimes of co- 
existence of bosons and fermions in the same space are 
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within the reach of experimental capabilities, and there 
is the exciting prospect of creating degenerate fermions 
in one dimension. Also we have shown that the low en- 
ergy density-density response of free fermions is identi- 
cal in form to that of weakly interacting bosons; this 
means that binary mixtures will have similar normal 
modes regardless of whether the components are bosons 
or fermions. These modes are analogous to zero sound 
modes and are essentially undamped at low momenta. 
Away from the low energy regime, the similarity of re- 



sponse functions does not hold, and the fermionic re- 
sponse in ID acquires interesting features, in particular 
a logarithmic divergence at twice the Fermi momentum 
which leads to periodic density modulation of the system 
analogous to Friedel oscillations. 
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